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Abstract 

If a catenoid is inverted in any interior point, a deflated compact geometry is obtained 
which touches at two points (its poles). The catenoid is a minimal surface and, as such, 
is an equilibrium shape of a symmetric fluid membrane. The conformal symmetry of the 
Hamiltonian implies that inverted minimal surfaces are also equilibrium shapes. However, 
they exhibit curvature singularities at their poles. These singularities are associated with 
external forces pulling the poles together. Unlike the catenoid which is free of stress, 
there will be stress within the inverted shapes. If the surface area is fixed, reducing the 
external force induces a transition from a discocyte to a cup-shaped stomatocyte. 

PACS: 87.16.Dg, 02.40. Hw 



1 Introduction 

Biological membranes are fluid bilayers whose mechanical properties are described on mesoscopic 
scales by a bending energy quadratic in extrinsic curvature [U El E]. If bilayer asymmetry or 
constraints are ignored, this energy is invariant under the global conformal transformations of 
three-dimensional space [HOE]. Thus the bending energy is not only independent of the size 
of the membrane, it is unchanged under inversion in any point, x — ► x/|x| 2 . As a consequence, 
every equilibrium configuration has a counterpart related to it by inversion which is also an 
equilibrium. Remarkably, this apparently unnatural symmetry does have physical consequences. 



In practice, what one does is pair inversions in special conformal transformations which can 
be treated perturbatively. The study of inversion in its own right, however, lies beyond the reach 
of perturbation theory. If the point of inversion lies on the surface, this point gets consigned to 
infinity; remote points on the surface, on the other hand get mapped into a neighborhood of 
the origin i.e. shapes do not simply get distorted, they may find themselves distorted beyond 
recognition; even their topology might change. Two shapes related by inversion will typically 
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describe very different physical conditions. Inversion has the potential to generate curvature 
singularities where none existed. 

In this paper, we will examine various consequences for the physics of fluid membranes of 
this non-perturbative aspect of conformal symmetry. To keep things simple we will look only 
at axially symmetric shapes. If our objective had been to identify compact isolated geometries 
without internal or external constraints, this symmetry would be a little too confining: the only 
equilibrium geometries of this kind are spheres and Clifford torii with ratio of the wheel to tube 
radius of y/2. 

The catenoid is also an equilibrium geometry albeit an infinite one. Finite sections of 
catenoid play an important role as a geometrical component of fluid membranes [lj. Because 
mean curvature vanishes, the contribution made by the catenoid to the bending energy is 
entirely topological, depending only on the boundary behavior. The central region of the 
catenoid provides the neck geometry permitting membranes to bud without incurring the large 
energy penalties one might expect to be associated with the high curvatures which form in the 
neck region. While the curvatures may be large, the mean curvature and thus the added local 
energy vanishes. This is also the geometry involved in the assembly of topologically complicated 
membranes. Thus it should not be surprising that such structures are abundant in biological 
cells. 

Spheres and torii are not very interesting under inversion. The catenoid, however, maps not 
only to a completely different equilibrium geometry, but to one that is also interesting physically. 
In particular, if the catenoid is inverted in its center of symmetry it transforms into a compact 
biconcave geometry: a discocyte. The two distal regions of the catenoid get mapped to the 
point of inversion. Thus the north and south poles of this geometry touch and, because the 
geometry possesses a tangent plane at these points, it is appropriate to think of the topology 
as spherical. 

The problem is: how does one interpret such an equilibrium? The only regular equilibrium 
geometries with spherical topology are round spheres. Without constraints, one would thus 
expect this geometry to equilibrate by inflating into a sphere. Physically, the only way such 
a deflated equilibrium could exist is if there are external forces counteracting the tendency to 
inflate. Such forces will act as a source of stress in the membrane [7j. The appropriate way to 
quantify this connection is to describe the geometry in terms of a conserved stress tensor [HI E] ■ 

Unlike a catenoid formed by a soap film spanning two rings which is under tension, a catenoid 
formed by an ideal fluid membrane is free of stress. However, the equilibrium described by its 
inverted counterpart is very different: it is under stress. As shown in [TU], the source of this 
stress is a pair of equal and opposite localized external forces at the poles pulling them together. 

External forces manifest themselves in curvature singularities. Despite the apparent smooth- 
ness of the inverted catenoid at the poles, there is a singularity lurking at these points. The 
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two principal curvatures diverge logarithmically as they are approached. 

What happens when the point of inversion is changed is also interesting. By translating 
this point along the axis of symmetry we break the up-down symmetry. The discocyte morphs 
continuously into a stomatocyte. It is also possible to describe this transition analytically. 
Surface areas, however, are not preserved under inversion. To follow the transition in a physically 
meaningful way, the geometries should therefore be rescaled so that they possess the same area. 
One may then track the volume as a function of the external force tethering the poles together. 
The maximum external force consistent with this constraint occurs for the symmetric discocyte; 
a higher value presumably would rupture the membrane. For each value below this maximum, 
there is a unique axially symmetric geometry. As the force is reduced, the discocyte first inflates 
asymmetrically until a point is reached where it is identifiable as a stomatocyte; thereafter it 
deflates into an inverted sphere within a sphere. We have not conducted a complete stability 
analysis; however, there are reasons to expect this description to be stable with respect to small 
perturbations breaking the axial symmetry. 

The paper is organized as follows: in section 2 we describe inverted minimal surfaces focusing 
on the stresses within them; in section 3 we provide a detailed physical interpretation of the 
inverted catenoid. Some of these results were announced in [10J. We finish with a brief 
discussion and an outline of our plans for future work. 



2 Willmore energy and inverted minimal surfaces 

A parametric description of an embedded two-dimensional surface in three-dimensional space 
is provided by the mapping (w 1 ,^ 2 ) — > X(m 1 ,m 2 ). The Willmore energy associated with this 
surface is given by 

1 „ fc „N /„ 1 



H[X] = J dA [K ab - -g ab K^j [K ab - -g ab K 



where g ab is the metric, K ab is the extrinsic curvature tensor, and K = g ab K ab is its trace 
(twice the mean curvature). dA is the induced element of area. This notation is summarized 
in Appendix I. H[X] is a measure of the energy associated with bending. It has the remarkable 
property that it is invariant with respect to the conformal transformations of the ambient 
space [4j: thus H\X] is invariant not only under the transformations of the surface induced 
by the Euclidean motions, translations and rotations, it is also invariant under transformations 
preserving angles: dilations X — > AX and inversion, 



Modulo a topological contribution proportional to the Gauss-Bonnet invariant H coincides 
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with the Helfrich energy 

H\X\= l -jdAK\ (3) 
There is only one local two-dimensional bending energy. 

The equilibrium surfaces of the energy satisfy the Euler-Lagrange equation 

- V 2 K + l - (Kg ab - 2K ab ) K ah K = , (4) 

where V 2 is the surface Laplacian compatible with the metric g ab . It is clear that minimal 
surfaces, satisfying K = are solutions of Eq.(|4]). However, because of the conformal symmetry 
of the shape equation, these surfaces map under inversion to new solutions of Eq.(|4]). This is 
true not just infinitesimally but also for finite conformal transformations. In particular, it is true 
for inversion. 

Let us first examine the behavior of minimal surfaces under inversion. In particular, let 
us first identify the equation satisfied by the transformed surface. It is well known that the 
principal curvatures transform under inversion by Cj — > Cj, where (see, for example, |15j ) 

<3 J = W (c,-2^). (5) 

Thus the transformed mean curvature K = C\ + C 2 is given by 

K = \X?(k-A<Z^). (6) 



X 



We conclude from Eq.© that a surface satisfying 



is mapped under inversion to a minimal surface, and conversely. It is clear that the only 
minimal surfaces mapping to minimal surfaces are the planes through the origin If the origin 
itself lies on the surface, K will diverge there unless the surface aligns along X sufficiently fast 
as |X| — > 0. While this is true if the original minimal surface has planar ends it is not true if the 
ends behave like a catenoid. Inverted minimal surfaces typically possess curvature singularities. 

Inverted minimal surfaces also satisfy Eq. (@]). Thus any solution to Eq.flT]) is also a solution 
to Eq.fllj). It is straightforward, but instructive, to confirm this fact with an explicit calculation: 
we note that, modulo Eq. (GO), we have 

4 ,„ 1 



d a K=—(K ab --g ab K) (X-e 



1 Eq. ((7|) can also be written as n- [V 2 X + 4X/|X| 2 ] = 0. Note, however, that there are no solutions consis- 
tent with peeling off n. Surfaces satisfying this equation can also be interpreted as the stationary configurations 
of the functional #i[X] = / dA/\X.\ 4 . 
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This makes use of the Weingarten equations, d a n = K ab e b . Eq.()H) follows (except possibly at 
the origin) on using the Gauss equations, V ' a e b = —K ab n as well as the contracted Codazzi- 
Mainardi equations, 

V a K ab - V b K = . (9) 

Here V a is the surface covariant derivative compatible with g ab . Unlike minimal surfaces, the 
equilibrium of the inverted counterpart involves a non-trivial cancelation between the Laplacian 
of K and the terms cubic in K ab . If the origin lies on the surface, however, equation ([7]) will not 
generally hold at this point and there will be a geometrical singularity there. In the following 
section we will show how to interpret this apparent pathology physically in terms of a localized 
external force acting on the membrane at this point. 

2.1 Conserved stress tensor and Noether charges for the inverted 
minimal surfaces 

The physical interpretation of geometrical singularities in the transformed geometry at the origin 
of inversion is facilitated by identifying a stress tensor with the geometry. This tensor is given 
in terms of the extrinsic curvature tensor by [SJ [5] 

P = K{K ab - ^g ab K)e b -d a Kn. (10) 

In equilibrium, f a is conserved so that V a f a = 0. It is simple to check that the local conservation 
law reproduces the shape equation (0J. 

The transformed stress tensor is given by 

P = |X| 4 (JX| 2 P + 4 ^K ab - ^9 ab K^j f 60 ) , (11) 

where 

f a = (e a -X)n- (n- X) e a . (12) 

The role of fg in the context of conformal symmetry has been noted elsewhere [16] . Curiously, 
it also provides an effective surface stress tensor for any external Laplace force [17J. 

In a minimal surface with K = the stress vanishes. However, it does not vanish in its 
inverted counterpart satisfying Eq. ([7]). The corresponding stress tensor is given by 

P = ( K<lb ~ \ gahK ) ffe0 ' (13) 

It is also straightforward to demonstrate that f a given by Eq. (fj~31) is conserved if X ^ 0: one 
uses Eq. ([5]) as well as the identity, 

Vhfoa = Gab**- + K b c ((X ■ e B ) e c - (a c)) , (14) 
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satisfied by ig. 

As noted in [10J, a curvature singularity at a point will be associated with a source of stress 
at that point. Let T be any closed contour on the surface of the membrane encircling the point. 
Stoke's theorem applied to the conservation law implies that the closed line integral 

£dsl a i a (15) 

is a constant vector F along contours that are homotopically equivalent to T on this surface 
[7]. Here 1 = l a e a is the normal to T tangent to the surface, and ds is the element of arc- 
length along T. The set of values of F are the Noether charges associated with the translation 
invariance of the energy. 

If the conservation law is valid everywhere, F must vanish on any topologically trivial loop. 
This is the case for a membrane with spherical topology. Thus, if F ^ on such a loop, there 
must be a source of stress within it [7J. In particular, a distributional source of stress manifest 
itself in a curvature singularity which is picked up by the line integral. F is a reparametrization 
invariant measure of the strength of singularity. 

Note also that, in equilibrium, the total external force acting on the membrane must vanish. 
In an axially symmetric membrane the forces operating at the poles must be equal and opposite. 

A concrete implementation of these ideas will be developed in the following section. 



3 Inversion of catenoids 



Conformal invariance permits a remarkably simple construction of a two-parameter family of 
axially symmetric solutions of Eq. (0}. Begin with a catenoid which solves the shape equation 
trivially as a minimal surface with K = 0. 

For a fixed axis of rotation, a catenoid £(.r ,£ ), is given as the level set $(i?, Z) — of the 
function 

$(R,Z) = R- — cosh^o^ + Co) ■ (16) 

It is characterized by a scale l/-Ro and an offset £,o/Ro along the axis. Thus any two catenoids 
are related by scaling and translation. They are the only non-trivial rotational minimal surfaces. 

The image of the catenoid S( Ro g \ under inversion in the origin is the surface S(^^ ) 
described as the level set 

wTz^ ~ -k cosh {wf&) + e °) = ' (17) 

2 Scaling the catenoid with an inverse length will give an inverted geometry with the 'correct' dimensions. 
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Figure 1: Section through symmetric biconcave discocyte 



This is a simple transcendental equation in the variables R and Z . 

One particular solution, the symmetric biconcave discocyte with £ = 0, is illustrated in 
Fig.(l). The geometric profiles for several different values of £ are displayed in Fig. (2). Sur- 
face area is not preserved under inversion; thus, as we discuss below, these profiles have been 
rescaled so that they all possess the same surface area. 

We note the following properties of these new surfaces: 

• The geometry is bounded for all values of £ - The circle R = 1/-R cosh£ on the plane 
Z = is mapped to the circle R = Rq/ cosh£ on the same plane. The remote regions 
of the catenoid map to the point of inversion. Thus the north and south poles of this 
geometry touch. As we will see, they do so with a common tangent plane Z = 0. In this 
sense, the inverted surface does not self-intersect. The topology is spherical. However, 
there will be a singularity at the origin. 

Generically, the behavior at the origin will not be so simple. Consider, for example, the 
fate of another simple geometry, the hyperboloid of revolution, under inversion in the 
origin. Asymptotically, the hyperboloid is conical. Because cones map to cones under 
inversion in the apex, the inverted geometry will be conical at the origin; it is not spherical. 
The existence of a tangent plane in the inverted catenoids is a consequence of the faster 
asymptotic growth rate (exponential) of the catenoid. 

• When £ = 0, the surface is a biconcave discocyte possessing up-down symmetry. This 
symmetry is a consequence of the symmetry of the original catenoid with respect to the 
point of inversion. 
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Figure 2: Geometric profiles for £o = 0, 1, 2, 3, 50, 100. The surface is generated by rotating the 
profile about the vertical axis. The area is normalized. 



• The symmetry is broken by translating the point of inversion along the axis. Beyond 
some critical value, £q = 1.86823 discussed below, the lower concave region becomes 
vanishingly small: there is a transition from the symmetric biconcave discocyte to an 
asymmetrical stomatocyte. 

• In the limit £ —> oo, the stomatocyte degenerates into a completely deflated geometry 
consisting of a sphere within a sphere connected by an infinitesimal neck, a catenoid. It 
is possible to describe the complete sequence from discocyte to limiting deflated stoma- 
tocyte analytically. 

It should be remarked that, technically, the lower concave region persists for all values of £ ; 
the ratio £o/Ro determines the position of the center along the axis; thus no matter how large 
a value we take for £o> there always will be a region of the catenoid where Z is negative; this 
region will map to negative Z in the inverted geometry and will contain a concave component. 
Of course, as the profiles above clearly illustrate, while this concave region is always present, it 
becomes small in comparison with the overall size of the geometry. How small will be described 
quantitatively in our discussion of the isoperimetric ratio. 

3.1 The isoperimetric ratio of the inverted shapes 

To examine the geometry quantitatively, it is useful to introduce a parametric representation 
of the catenoid in terms of the angle that the tangent to the meridian makes with a plane of 
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constant Z: 

j R(9) = ^csce, Z(S) = In tan 0/2 + . (18) 
Rq Rq Ro 

where lies in the interval [0, 7r]. This angle also provides a parametrization of the inverted 
catenoid E: 

5 ro % RpsmQ - # o sin 2 0(lntanf + £ ) 

1 J l + sin 2 0(lntan§ + £ o ) 2 ' [ ' 1 + sin 2 (In tan f + £o) 2 ' 1 ' 

The interval of does not change under inversion. A word of caution: and its inverted 
counterpart are not the same angle i.e. has nothing to do with the angles which are 
preserved under conformal transformation H. While it might be more natural to parametrize the 
inverted surface by 0, the advantage has to be weighed against the complicated form of the 
corresponding functional forms of R and Z. 

A straightforward calculation shows that the parametrization by is isothermal. The metric 
on the inverted catenoid is described by the line element 

ds 2 = R 2 Q n 2 (0) [d& 2 + sin 2 dip 2 ] , (20) 

where ip is the polar azimuthal angle, and the conformal factor Vt is given by 

/ \ 2 

^(0)^ = 1 + sin 2 log tan- + £ . (21) 



The area of the surface is given with respect to this parametrization by A = AuRqA^q) where 

A (&) = 1 r dB sin Q (0) 2 , (22) 
2 Jo 

with range (0, oo). The area of the surface is used to express R as a function of £ for 
a fixed area A as R (£ ) = \J A/A-nA (£o) The volume enclosed by the surface is given by 
V = ^ V(f ) , where 

V{£a) = \[dQ sin 3 e(l + cos© (log tan®+fo)) fi 3 (0) , (23) 

and thus the isoperimetric ratio v, defined by v = V/Va = 3y/4arV/A^ where Va is the volume 
of a sphere of surface area A, is given by 

^o) = V(£ )A4(£o) § • (24) 

3 More generally, e a ■ k ^ e a • k, where k is a fixed unit vector. 
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The maximum value, v = 1 is attained by a sphere. Clearly the constraint that the poles touch 
introduces bending energy. As a consequence the value v = 1 is not attained in any equilibrium 
geometry with touching poles. 

The isoperimetric ratio is plotted as a function of £ in figure 3(a). v (£ ) possesses three 
extrema, two maxima and one minimum. The maxima t> Max = 0.66564 occur when Q = 
±1.86823. The corresponding geometry marks the transition from discocyte to stomatocyte 
(depicted in figure 2). The shallow local minimum t> M m = 0.647639 occurs when £ = 0, 
the symmetric biconcave geometry. We observe, in particular, the interesting fact that the 
maximally inflated geometry is not the symmetric discocyte. v vanishes asymptotically in the 
limits £ — > ±oo. This property is evident in the corresponding profile (see figure 2). 

We observe, in particular, that i>(£o) may not be inverted for £ on the interval £ G [0, oo). 
This indicates that it cannot be used in place of £ as the 'order parameter' to describe the 
discocyte-stomatocyte transition. As we will see, the appropriate physical parameter is the 
Noether charge or external force F. This physical interpretation will be developed in the 
following sections. 



Reduced Volume Area Ratio 




5 10 15 20 -10 -5 5 10 

ft ft 



Figure 3: (a) Isoperimetric ratio vs (b) Area ratio vs £o- 



It was pointed out in the previous section that, from a technical point of view, the geometry 
remains biconcave for all values of £ . We will now demonstrate that one of the concave regions 
becomes vanishingly small when £ is larger than the critical value Q = 1.86832. To do this, 
we determine the fraction of surface area lying above Z = 0. We thus define the ratio 



A 



r d& sinefi(e) , (25) 
Je (Zo) 



where 6 (£o) = 2 arctan (exp (— £o)) is the tangent angle along the meridian where the surface 
intersects Z — 0. For example, when £ = 0, 6* (0) = f and the area ratio is one half. Figure 
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3 (b) illustrates clearly the rapid approach of the ratio to 1 above the critical value, Q. 



3.2 Curvature singularities 

We expect a curvature singularity at the poles. This is because the only regular equilibrium 
geometries with spherical topology are the round spheres. We will now confirm that the ge- 
ometries described by Eq. ( TlTl) are regular everywhere except at the poles where they display 
a (logarithmic) curvature singularity. The strength of this singularity will be related to the 
external force tethering the poles together. 

In the case of the symmetric discocyte it is very simple to study this singular behavior 
analytically. To do this, note that, in the neighborhood of the origin, Eq. ( TTTT ) is approximated 
by 

ex P -52" ( 26 ) 



R 2Rq V R? 
which can be inverted for Z as a function of R 



Z ~ ^V£) (27) 



Ro \RqJ \2Rq, 

when R « 0. Whereas Z tR vanishes at R = 0, it is clear that both the curvature along the 
meridian C±_ ~ Z RR , as well as that along the parallel, Cm ~ Z R jR, diverge logarithmically 
(see appendix II). Indeed, C\\ and C± exhibit identical logarithmic divergences: 

C h C ± « -^-log^ H 



Ro V 2R , 

The poles are umbilical points of the geometry, albeit in a singular way. 

We will now show that this qualitative behavior holds for all values of £ - However, a different 
strategy becomes more appropriate when £ ^ 0. We have the following exact expressions for 
the principal curvatures with respect to the parametrization given by (Tfj 

dlR dp,Z - dlZ dp,R 



3 

\2\ 2 



((d @ R) 2 + (d @ Zf 

— f3 + 2cos6 (^logtan® + ^o) + sin 2 9 (^log tan ® + ^o) J, (2f 



and 

Ch 



9« Z 



R((d e R) 2 + (d e zy 



2\ 2 



~ (l + 2 cos (log tan | + ^o) -sin 2 9 (log tan | + ^o) ^ • (29) 
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It is clear, by inspection, that the two principal curvatures diverge logarithmically at the poles 
6 = and 6 = it. Note, however, that their difference remains finite with the value C± — C\\ = 
—2/Rq as consistency with Eq.(2) of [10J demands. The angular dependence of the curvatures 
for different values of £ is plotted in Figures 4 and 5. 




Figure 4: Ci vs. G for three different values of £ : Ro = 3, £ = (left); R = 3^(0) /A{2), 
£o = 2 (center); R = 3^/A{0) /A{3), £ = 3 (right). 




Figure 5: C\\ vs. G: R = 3, £ = (left); R = 3^-4(0) /A{2), £ = 2 (center); R = 
3^A(0)/A(3), Co = 3 (right). 



Note, in particular, how sharp the divergence is in the neighborhood of the south pole. Note 
also that both curvatures diverge to infinity. This peculiarity is associated with the logarithm; 
it is not apparent visually in the corresponding profiles due to the existence of a tangent plane 
(figure 2). 

The singularity at the poles indicates that these surfaces, unlike the catenoid, require external 
forces along the axis of symmetry to support them. 

It should be pointed out that the characteristic logarithmic singularity we have examined 
is, in fact, predicted by the linear theory. Consider the linear approximation to the Monge 
representation of the surface in terms of its height h above a plane. If the gradients of h are 
small, the bending energy H can be expressed as 

H= l -fd 2 x(Vlh) 2 + 0(A) (30) 
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where Vj_ is the gradient on the plane. The corresponding Euler-Lagrange equation is the 
biharmonic equation equation on the plane, — (V^) 2 /i = 0. It is straightforward to show that 
the most general axially symmetric solution of this equation is given by 

h = c + c x R 2 log R 2 + c 2 log R + c 3 R 2 (31) 

The linear Monge approximation to the inverted catenoids is given by solutions with c ,c 2 = 0. 
Of course, it is beyond the scope of this approximation to model their global behavior. 



3.3 Singularities as manifestations of sources of stress 

Axial symmetry dictates that the force on the membrane must be directed parallel to the axis, 
F = —2irck. We have shown recently that c is the constant appearing in Eq. (2) of [TU] , 
Thus using this equation and Eq. ([7]) we find that an inversion of the catenoid must satisfy the 
following equation 

Cx-C,^. (32) 

If we now use the explicit expressions for the principal curvatures we identify c = —A/R . This 
is a reparametrization invariant measure of the singularity. It is also possible to see this directly 
using the expression given in Eq. ( Tl5l) : without loss of generality, choose the contour to be a 
closed circle of polar radius R encircling the north pole. Now 1 is the tangent to the meridian 
and ds = Rd<fi. We do not use Eq. ffTBl which is ill-defined at the pole. Using ( lTg|) and 
and that on the north pole n = — k, we obtain 

lj a = i(C|-Cf)l-(C ± + q|)'n 



{C ± + C\\)'k 



(33) 



As a consequence, 

F = 8vr.Ro k , (34) 
so that c = —A/R . The sign indicates that the force is always directed towards the interior. 

3.4 The Phases of a tethered membrane 

We now possess a physical interpretation of the inverted geometry as an equilibrium with a 
pair of equal and opposite localized external forces applied at the poles holding them together. 
This might be compression applied outside, or the tension supplied by an internal tether. The 
curvature singularity at the poles is a manifestation of local forces. It is unrelated to the 
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fact that they touch: they touch simply because all distant parts of the original catenoid get 
mapped to the origin, a peculiarity of the exact solution we have been handed. In a subsequent 
publication we will construct equilibrium configurations with a fixed non-vanishing interpolar 
distance. 

The axial force c may be interpreted as the order parameter describing the transition from 
discocytes to stomatocytes. The axial force and the reduced volume may be used to characterize 
axially symmetric compact shapes. The geometries with tethered poles that we have considered 
define a finite trajectory on this space illustrated in figure 6. Each point on this trajectory 
represents a specific geometry: the discocyte lies at one end (with the maximum fore c\ = 
—4/3), the limiting geometry consisting of a sphere within a sphere at the other (vanishing 
force c = 0). The maximally inflated equilibrium geometry (c* = —0.711923) provides a 
natural division of the trajectory into two phases: discocytes and stomatocytes. 



4 Discussion 

The invariance of the bending energy of two dimensional surfaces with respect to inversion in 
any point has been used to obtain exact results for a family of deflated biconcave equilibrium 
geometries with spherical topology and with tethered poles. The external force tethering the 
poles is imprinted as a curvature singularity on the membrane geometry. We end up with a 
geometry describing a completely different physical setup from the catenoid geometry we began 
with. 

There is a, of course, a superficial similarly to the discocyte-stomatocyte transition of fluid 
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membranes induced by a change in the bilayer asymmetry. The quantitative study of the 
discocyte shape of a deflated fluid membrane and its stability has a long history dating back 
to Deuling and Helfrich's pioneering work in the seventies [llj followed, after a long pause, 
by a series of ground-breaking studies by several groups (see, for example, [12] as well as 
[21 [T3] [14] for reviews). The transition explored in this work is associated with a change in the 
bilayer asymmetry. An energy penalty or constraint on the area difference is introduced. These 
constraints break the conformal symmetry of the bending energy pointwise. In the inverted 
catenoid, on the other hand, there is no bilayer asymmetry. The transition is associated with a 
change in the localized external forces. The action of these forces breaks conformal symmetry 
only at isolated points. As a consequence, it is possible to constrain the isoperimetric ratio 
without sacrificing the conformal invariance of the energy. While it is tempting to emphasize 
the analogy between the two systems, it is probably a mistake to dwell on it. 

These differences notwithstanding, the simple model we have described is relevant to the 
study of biological membranes in other ways. For the action of forces external to the membrane 
controlling its shape is the rule rather than the exception in biology. The cytoskeleton pushes 
and pulls on the plasma membrane [IB]. Indeed, our toy model might after all be relevant to the 
stability of the biologically relevant biconcave discocyte shape assumed by a healthy erythrocite 
which involves this interaction in an essential way. 

It is also possible that our model is relevant for the theoretical study of the geometry of 
two very important organelles: the Golgi complex and the mitochondrion. It is known that the 
Golgi complex is not an isolated equilibrium structure [!§]. It is stabilized against breakup by a 
complex network of microtubules that is dismantled and reassembled during mitosis. And the 
convoluted inner membrane of the mitochondrion is likely to be under compressive stress due 
to its confinement within the outer membrane [20J. 

In this context it is perhaps worth stressing that, in principle, it is always possible to recon- 
struct the distribution of external forces {Fi, . . . , Fjv} from a knowledge of the geometry in 
the neighborhood of an appropriate set of contours. The external forces leave their imprint on 
the local membrane geometry. As imaging techniques improve, this could become a very useful 
diagnostic tool in the study of intracellular membranes. 

The solutions we have described are very special in that the poles touch. Intuitively, one 
would expect that solutions exist with poles held a finite distance apart. In the early nineties 
Ou-Yang Yhong-can, Naito and Okuda identified an exact analogue of such configurations for 
a model with a finite spontaneous curvature [22J. Generally, however, the effect of an external 
force will be to pull out two cylindrical tethers. Tether pulling in the context of a physically more 
realistic Hamiltonian has been the subject of some very nice recent theoretical analysis (see, for 
example, [23, 24, 25]). Such geometries will occur as solutions of the axially symmetric shape 
equation we derived in reference [10] . A detailed description of the solutions of this equation 
is, however, beyond the scope of this paper and will be addressed in a forthcoming paper. It 
will come as no surprise that conformal symmetry has a role to play in their characterization. 
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We have said nothing about the stability of the shapes we have discussed. Expanding the 
bending energy out to second order in deformations about a minimal surface, one can show 
that [26J 



where $ = <5X ■ n is the normal deformation and 1Z is the scalar curvature (twice the Gaussian 
curvature CiC 2 ). Remarkably, this expression depends only on intrinsic geometry. It should be 
compared with the corresponding expression for the area, 



In general, for any minimal surface 1Z < 0. Thus the catenoid spanning two rings becomes 
unstable as a minimal surface beyond a certain maximum separation. By contrast, note that 
the expression given by ( 1351) is manifestly positive. Thus, if tension may be ignored, the same 
catenoid is stable as a Willmore surface. The conformal invariance of the energy should hold 
order by order in perturbations theory. This suggests that the inverted shapes are also stable. 
While one does need to confirm that there is no subtlety associated with singularities, conformal 
invariance clearly simplifies the analysis of stability. 

So much for axially symmetric configurations. It is clearly possible to generate non-axially 
symmetric geometries held by pair of tethers by inverting these geometries in any point off the 
axis of symmetry. It appears reasonable to conjecture that all equilibrium shapes with a pair of 
tethers are generated by inversion. 

It is also clear, however, that there is much more to this story. For the catenoid is but the 
simplest of a vast and growing number of interesting minimal surfaces. How do we describe 
equilibria held in place by three or more tethers? It is not unreasonable to guess that the 
inversion of non-axially symmetric minimal surfaces will play a role. We are currently examining 
the properties under inversion of the natural generalization of the catenoid, the n-noids of 
Jorge and Meeks [27]. This could provide a non-perturbative handle on yet another problem of 
interest in membrane biophysics: how does one identify the stable configurations assumed by 
proteins when the interaction between them is mediated by membrane curvature [28J? 

Appendix I 

The tangent vectors adapted to the parametrization are e a = <9 a X and n is the unit normal. 
The induced metric on the surface and the extrinsic curvature are given by g ab = e a ■ e& 
and K ab = e a • c^n respectively [29] ETJ- Indices are raised with the inverse metric g ab . 
The element of area is given by dA = \/det g ab d 2 u. K denotes (twice) the mean curvature: 




(35) 




(36) 



K = g ab K ab 



Appendix II 
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The principal curvatures are given by 

C ± = G', C\\=nnG/R 

where denotes the angle that the tangent along the meridian makes with the planes of 
constant Z, and the prime denotes a derivative with respect to arc-length along the meridian. 
The arc-length is given by 

ds 2 = (1 + Z 2 R )dR 2 

Thus cos© = R' = -1/(1 + Z 2 R fl 2 « -1, sin = -Z' = Z :R /(1 + Z 2 R fl 2 « Z >R , and 
©' « Z tRR . 
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